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fT^ ' Abstract The occurrence of aftershocks following a major financial crash mani- 

fests the critical dynamical response of financial markets. Aftershocks put additional 
stress on markets, with conceivable dramatic consequences. Such a phenomenon 
L^ ' , has been shown to be common to most financial assets, both at high and low fre- 

ry\ ' quency. Its present-day description relies on an empirical characterization proposed 

by Omori at the end of 1800 for seismic earthquakes. We point out the limited pre- 
dictive power in this phenomenological approach and present a stochastic model, 
based on the scaling symmetry of financial assets, which is potentially capable to 
O^' predict aftershocks occurrence, given the main shock magnitude. Comparisons with 

S&P high-frequency data confirm this predictive potential. 
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1 Introduction 

It is not uncommon for financial indexes or asset prices to experience exceptionally 
large negative or positive returns which trigger periods of high volatility, the case 
of abnormal negative returns corresponding to market crashes. An understanding of 
the dynamical response of the market to a main shock is of great interest because it 
may help, e.g., in the definition of emergency plans for financial crises, or for risk 
management. 

There is a clear analogy between the behavior of volatility after a main financial 
shock and that of the seismic activity after an earthquake of exceptional magnitude 
in geophysics lfT2l . Omori ifTOl , with a subsequent modification by Utsu |fT6l , estab- 
lished an important empirical law describing the frequency of occurrence of seismic 
events above a given threshold after a main earthquake. The characterizing feature 
of this law is the decay as a power of time, t, of the rate of occurrence of after- 
shocks above the threshold, indicating the absence of a characteristic time scale in 
the manifestly non-stationary Omori regime. More precisely, according to Omori 
the number, n(f ), of aftershocks per unit time above a given threshold Ga is given by 

n{t)=Kit + T)-P, (1) 

where K, T,p depend on the aftershock threshold Ga, and also on the specific mag- 
nitude of the main shock earthquake. Equivalently, the Omori law can be expressed 
in an integral form as 

Nit)^-^[it+Ty-p-T'-p] (2) 

1 — p 

if p 7^ 1, or N{t) = K ln(f/T +1) if p = 1, where N{t) is the cumulative number 
of aftershocks up to time t after the main shock. Lillo and Mantegna Ul were the 
first to verify the validity of an analog of the Omori law for the volatility in Finance 
after a main crash. They also showed [8] that standard dynamical models of index 
evolution, like GARCH, are not adequate to reproduce financial Omori-like regimes. 
Several studies ||7] |8] [T3] [141 [HI HI [HI verified the presence of Omori regimes 
under various market conditions, triggered by financial crashes li7l[8l [T3l[T4l[T8l . by 
volatility shocks Q, and even by U.S. Federal Open Market Committee meetings 
lUn . In particular, the Omori law in finance has been upgraded to a more general 
characterization of market dynamics by Weber et al. iTTSl , who pointed out that this 
law holds on a wide range of time scales, with aftercrashes of a main shock playing 
the role of main crashes for even smaller aftercrashes, etc. 

The above mentioned studies make clear the connection between financial Omori 
processes and long-range dependence of the volatility. They also show that a mod- 
ulating, time dependent scale for the returns must be considered in order to account 
for the manifest non-stationarity of the Omori process. At the same time, they em- 
phasize the limits in the predictive value of the Omori law. For example, the pa- 
rameters K and T need to be adjusted for each aftershock threshold considered (See 
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below). This holds also for the exponent p of the power law decay, which should 
be expected to be the most robust parameter In addition, there is no idea of how 
the parameters could be linked to the magnitude of the main shock. These limits 
reflect a lack of adequate modeling for the dynamics of financial indexes, especially 
in regimes like those covered by the Omori law. In recent contributions ElIU, some 
of the present authors have proposed a model for the dynamics at high frequency 
of exchange rates or stock market indexes, which takes into account most of the 
relevant stylized facts. Among them, the martingale character of index evolution, 
the manifest non-stationarity of volatility detected in well defined daily windows of 
trading activity, the anomalous scaling properties of the aggregate return probability 
density function (PDF) in the same windows, and the strong time autocorrelation of 
the elementary absolute return. This model for high-frequency data, which applies 
more general ideas about the time evolution of financial indexes 121 [15] IS, has also 
been tested JS) by comparing its predictions with the statistics of ensembles of daily 
histories all supposed to reproduce the same underlying stochastic process. It has 
been also shown lH that some arbitrage opportunities revealed by the model could 
be successfully exploited by appropriate trading strategies. 

In the present contribution we address the problem of describing with such a 
model the Omori processes which may be detected within these daily windows. Our 
goal is to show that, after proper calibration, this model allows the prediction of the 
aftershock rate within an Omori regime, given the value of the main shock magni- 
tude. Indeed, we provide analytical expressions for the rate of financial aftershocks 
with explicit dependence on the magnitude of the main shock and on the aftershocks 
threshold. By comparing our predictions with high frequency data from the S&P 500 
index we show that these quantities are sufficient to determine the Omori response 
without further fitting parameters. Our success is partly due to the fact that we are 
able to identify the Omori processes within a context for which non-stationarity is 
well established [5 J and amenable to modeling 13] H) . In an interday context, the 
question of the applicabihty of the models of Ref. ID to Omori regimes has already 
been raised in Ref. Q. 

This note is organized as follows. In the next Section we briefly recall the model 
of Refs. [3 4j and present the procedure of calibration. In the third Section we 
discuss the selection of Omori-like processes from our database and show how our 
model can be used to analytically describe these processes. In the fourth Section we 
compare the results of the properly calibrated model with the statistical records at 
our disposal for the S&P 500 index. The last Section is devoted to general discussion 
and conclusions. 



2 Model calibration 

Let us consider the successive (log-)returns over ten minutes intervals of the S&P 
500 index S{t) for daily windows from 9.40 a.m., Chicago time, to 1.00 p.m.: 
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R,=lnSit+l)- lnS{t), f = 0, 1 , . . . , 19, 



(3) 



where the time is measured in ten-minute units and we have set f = at 9.40 a.mj^. 
A statistics made over the ensemble of 6283 available daily histories from 1985 to 
2010 shows PI that a stochastic process supposed to generate the successive returns 
Rt in a generic history of the ensemble is consistent with the following joint PDF: 



PRo.R,,...,R,i''0^''i,...,r,)^ / dap{a)Y[ 



exp 



'2a^ 



Ina^o 



where 

with D > 0, / = 0, 1, . . . , 19, and p (a) > with 



i: 



da p{a) = 1. 



(4) 



(5) 



(6) 



This joint PDF is a convex combination of products of Gaussian PDF's for each in- 
dividual return. The PDF p (a) weights this combination and introduces a nontrivial 
dependence of the returns from the preceding ones. For D ^ 1/2, the coefficients a, 
make the process increments non-stationary, and modulated by the exponent D. 

The calibration of the model can be done by direct comparison of its predictions 
with the main features of the PDF's of the 10-minute returns /?,'s, or, alternatively, 
with those of the aggregate returns Y!i=o^i ElSl. Here we follow the second option. 
Since the model predicts for the PDF of the aggregate return L/=o^i satisfaction of 
an anomalous scaling of the form 



pnw^My 



1 



where the scaUng function g is expressed as 



(7) 



?(r)= / dapia) — ^== , 



(8) 



one can determine D through a fitting of the power law f -dependence of the moments 
of pyt g . Indeed, for q E R, according to Eq. (|7]i 
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'n 


_ !=0 





= E[|/?orv 



^i fr(?) 



(9) 



' In order to keep contact with ordinary notations for the Omori law, in this paper we change 
slightly our usual conventions by shifting the origin of time by one unit with respect to, e.g., Refs. 

mm. 
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Fig. 1 Calibration of the scaling exponent D. The empirical values for 7(17) are reported using Eq. 
J9j as an ansatz (points). A linear regression for < ^ < 2 gives /(ij) = q D with Z) ~ 0.35 (dashed 
line). 



with Y{q) = q D, and provided that the moment E [l^oT'] exists. In Fig. [T] we report 
the empirical values for Y{q), using Eq. (|9]l as an ansatz. To calibrate D, we make a 
linear data regression for q <2, since for higher moments a multiscaling behavior 
iflTl is detected (See Fig.[Ili. The result is D ~ 0.35. 

A particularly simple expression for the joint PDF pRf,.Ri r, is achieved if the 

integration on a can be worked out explicitly in Eq. (|4]i. This is indeed the case if 
we choose an inverse-gamma distribution for a^ ^. Equivalently, we may set 



Pi<y) = 



2i-fj3' 



^(f) 



70:+ 1 



exp 



2C72 



(10) 



where the exponent a determines the long-range behavior of g according to g{r) ^ 
l/r"+' for \r\ ^ 1, and j3 is a scale parameter determining the distribution width. 
Performing the integration on a in Eq. (|4]i we obtain a multi-variate Student PDF: 
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0.02 0.04 



Fig. 2 Calibration of the parameters a and j8. The empiiical PDF's for Y,'i=o^i ^^ various t are 
rescaled according to Eq. ([T} with the previously calibrated D = 0.35 (points). The parameters 
a and j8 are then fitted using Eq. l IlU with / = (dashed line), yielding the values a = 3.5 and 
j8 = 2.9 • 10^-'. An upper bound to the empirical analysis is posed at the (J,„i,x = 0.02 for the 10- 
minute volatility. 



PRa,Ri....,R,iro,n,---,rt) 



j3«r 



' a+t+l ' 



n'-^r(^] 



i8^ + 4 + 4 



(11) 

As we will show in the following, an explicit form for Pr^.r^ r, enables us to obtain 
a simple analytic expression for N{t). Unlike in previous papers ||3] |4], we thus 
choose here the functional form in Eq. (fTOt for p . Besides D, the other parameters 
of the model, a and j5, are calibrated by first data-collapsing the empirical PDF's 
for Y!i=o^i according to Eq. (|7]i with D = 0.35, and then by fitting a and j5 on this 
data-collapse using Eq. (fTTT i with f = 0. The result is given in Fig. |2] In summary, 
the result of the calibration procedure is the triple (a,j3,D) = (3.5,2.9- 10^^,0.35). 
The ensemble of histories at our disposal is relatively poor. This implies, as can 
be appreciated in Fig. |2] that some rare events fall significantly out of the scaling 
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function, since a much larger number of histories would be needed to correctly char- 
acterize their frequency of occurrence. The multiscaling behavior shown in Fig. [T] 
could be at least partly related to this effect. The Omori events are precisely related 
to extreme events. In order to obtain a reliable statistics of the aftershocks, we im- 
pose thus an upper bound a„,ax to the absolute value of the returns 7?,'s included in 
our empirical analysis (See Fig.|2}. Once done this, the overall agreement of the em- 
pirical data with the various model predictions gives a convincing validation of the 
model itself (See also JU). Still, the agreement shown in what follows with respect 
to the Omori processes must be intended as a first important result, which calls for 
more extensive analysis also in terms of the calibration procedure. 



3 Aftershock prediction 

As already mentioned above, in the present analysis we are going to identify and 
select Omori processes, which are manifestations of non-stationarity, within a pro- 
cess which manifestly turns out to be with non-stationary returns in its ensem- 
ble of daily realizations. This is a simplification which marks an important differ- 
ence with respect to the problem of modeling the Omori regimes revealed in Refs. 
|I7][8][T3][14][18]|9][TT|, where they were extracted from single time series expected 
to be globally stationary on long time scales. In the perspective of our approach 
here, dealing with a process which is by itself time-inhomogeneous offers the ad- 
vantage that the selection of Omori processes does not imply the need of identifying 
how their non-stationarity emerges from an otherwise stationary global behavior 
In a version of our model suited for describing single, long time series of returns 
lfT5l[T9l , the necessity to consider random exogenous factors influencing the market, 
leads us to switch-on at random times some time-inhomogeneities formally similar 
to those characterizing the model of the previous Section. This is achieved by set- 
ting a; — I concomitantly with these random events (See also 12113] |6][T]). In such 
a context it is not a priori clear whether or not the start of an Omori process should 
imply putting a, = 1 in correspondence with the time t of the main shock. This dif- 
ficulty is also accompanied by the need of implementation of a more complicated 
calibration procedure llT9l with respect to the one presented here. 
As shown below, remarkable results of our analysis in this note are: 

1. that the selected processes are legitimately classified as Omori-like in the sense 
that they can all be fitted by the Omori law; 

2. that the description one obtains for them based on the model presented in the 
previous Section contains explicit dependencies on the intensities of the main 
shock and on the aftershocks thresholds. 

This endows our approach to the Omori regimes of a predictive potential which, if 
confirmed by further analysis, could be exploited by decision-makers under crisis 
conditions. 
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We select as Omori processes all those histories in the above S&P 500 en- 
semble for which the initial absolute return, |ro|, besides being smaller than Gmax, 
also exceeds a main shock threshold c,,,. At variance with the analysis in Refs. 
|l7l[8l[T3][l4l[T8]|9l[TT|, we consider, in place of a single time series, groups of his- 
tories for which (7„, < | ro | < Gmax ■ As far as the aftershocks are concerned, we record 
for each of these histories the elementary returns which exceed in absolute value an 
aftershock threshold (7„ and are below the main shock value |ro|: (To < |r,| < |ro|, 
for / > 1 . The parameter Oa is an important one to be fixed in any analysis of the 
Omori law. Again, by imposing the aftershock magnitude to be smaller than that of 
the main shock we reduce the influence of extreme events in our limited dataset. We 
decided to search for main shocks occurring right at the beginning of the daily time 
window described by our model for two main reasons. In first place the ensemble 
average volatility on 10 minutes intervals is maximal in the first interval. Secondly, 
a main shock occurring right at the beginning of the time window leaves the max- 
imum possible time for the development of the subsequent Omori process. While 
we will limit ourselves below to discuss such optimal case, different choices are of 
course possible. 

According to the above selection procedure of the Omori processes, the cumu- 
lative number of aftershocks N\,-q\ (t) after a main shock of magnitude |ro| is given 
by 



^k-o|(0 = E 



]Ma„<\R,\<\Ro\) I l^ol = ko| 



= 1 



(12) 



where 1(ct„<|R;|<|/;o|) is the indicator function, yielding 1 if (7„ < \Ri\ < \Rq\ and zero 
otherwise. Using Eq. (fTTT i. through a change of variable it is straightforward to show 



,=1 ^o„ PRoiro) 
'a+2\ t n ^ 



2 ^{-^ 



^^(^)-Wi^ 



^Z-'V^^rfxCl+x^)-^. (13) 



^'0 

If in the considered ensemble of histories there are M realizations 



L J m=l,2 M 

in which we register a main shock, i.e., a,„ < |^o" | < Omax, then the cumulative 
number of aftershocks A^(f ) is obtained through the sample average 

2 M 
^ m=l I" I 

where we stress the fact that each Ni (m)! (?) is conditioned to the main shock magni- 



'0 



tude 



J") 



. Notice that since with the available dataset the selected main shocks con- 
stitute a small sample (see next Section), we use here the sample average rather than 
the ensemble one to get the number of aftershocks conditioned to C7„, < \Ro\ < G,„ax- 
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Fig. 3 Fitting of the empirical aftershock at different thresholds o;, (points with error bars), with 
the Omori law in Eq. |2) (dashed lines). Fitted parameters are reported in Table [T] 

Table 1 Omori parameters in Eq. |2]fitted from the empirical data. 
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4 Comparison of the model predictions with the statistics of 
aftershocks 



Our choice for the thresholds <7m and (Jmax is such that the absolute first returns 
for which C7„, < |rg '| < (7„„„ are quite exceptional. They occur with 27/6283 ~ 
0.4% frequency in our ensemble; Only 3 realizations have |ro| > Omax and are thus 
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Fig. 4 Comparison between the analytical model predictions for different aftershock thi'esholds 
(Tj, (dashed lines) with the same empirical S&P data reported in Fig.[3](points with error bars). 



excluded. Accordingly, we analyze the averaged N{t) of aftershocks for these M = 
27 main shocks. A first point to clarify is whether the recorded rates are well fitted by 
the Omori law in Eq. (|2]i. This is shown in Fig.|3] where many sets of data for A^(r), 
obtained with different aftershocks thresholds (7„, are indeed fitted by the Omori Eq. 
(|2|i. In Table[T]one also realizes that K, T and p need to be varied for each (7„ in order 
to reach a satisfactory fit. In particular, by changing K and T it is even possible to 
obtain reasonable fittings also with p > I (See also |l9l)- This parameters variability 
makes it very difficult to use the Omori law to predict the aftershock occurrence for 
a given main shock magnitude and aftershock threshold. 

Model predictions on the same set of data fitted in Fig |3] are instead given in 
Fig. m Dashed lines in Fig. |4] are obtained on the basis of Eqs. ( 113114b with the 
parameters {a,j5,D) resulting from the calibration discussed in Section|2] The only 
difference among the curves is the value of the aftershock threshold (7„. The agree- 
ment of the analytical predictions with the data and the sensitivity of the curves 
to the variation of the aftershock threshold are remarkable. This shows that our 
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model potentially provides a satisfactory and parameter-free description of Omori 
processes. 



5 Conclusions 

We have shown here in the case of the S&P 500 index, that a model suited for the 
description of the high frequency market dynamics allows also to predict Omori 
regimes following exceptional extreme events. Within the class of events consid- 
ered, the model specifies the dependence on the main shocks intensities and on the 
aftershocks threshold. As such, its description goes far beyond the limits of the 
Omori phenomenological law. 

Besides providing a further validation of the model of Refs. |[3] ID, the results 
presented here encourage to extend similar analysis to cases in which the Omori 
processes are to be selected within a process which is globally stationary. For the 
modeling of these processes, our recipe 1T91 is that of switching-on at random some 
non-stationarities ascribable to coefficients like the a, defined above. Global station- 
arity of the process on long time scales is then guaranteed by the fact that empirical 
averages are in this case made by considering time intervals sliding along the sin- 
gle long history |fT9l . While it is conceivable that in many cases main shocks are 
localized close to resets of the time inhomogeneity, this is not true in general. Some 
attempts to strictly identify main shocks with restarts of them inhomogeneity in the 
model (fl, = 1) already gave some preliminary agreement with the data. A more 
general discussion is however needed [[T9l . 
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